Supersymmetric formulation of multiplicative white— noise stochastic processes 
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We present a supersymmetric formulation of Markov processes, represented by a family of 
Langevin equations with multiplicative white-noise. The hidden symmetry encodes equilibrium 
properties such as fluctuation-dissipation relations. The formulation does not depend on the par- 
ticular prescription to define the Wiener integral. In this way, different equilibrium distributions, 
reached at long times for each prescription, can be formally treated on the same footing. 
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I. LANGEVIN REPRESENTATION OF A 
MULTIPLICATIVE WHITE-NOISE PROCESS 



Several years ago it was recognized that some stochas- 
tic processes have a hidden symmetry [J-Q called super- 
symmetry (SUSY). From a mathematical point of view, 
SUSY shows up as an invariance under linear transforma- 
tions that mix commutative as well as auxiliary anticom- 
mutative variables in a path integral representation of the 
process. Physically, it encodes equilibrium properties of 
the system. Some of the constraints it imposes on correla- 
tion functions (Ward-Takahashi identities) are related to 
fluctuation-dissipation theorems [4]. This property has 
acquired a renewed interest due to the growing impor- 
tance of stochastic out-of-equilibrium systems 5] . In this 
sense, it is possible to understand out-of-cquilibrium dy- 
namics as a symmetry-breaking mechanism. SUSY prop- 
erties have been studied extensively for additive stochas- 
tic processes [1,0] and, more recently, for non-Markovian 
multiplicative processes [8]. However, the important case 
of Markovian multiplicative white-noise systems has re- 
mained elusive [9|. The main difficulty in the SUSY for- 
mulation of a Markov process resides in the great vari- 
ety of prescriptions to define the Wiener integral, which 
produces several stochastic evolutions with different final 
steady states. Moreover, time reversal transformations 
mix different prescriptions. In this letter we present a 
supersymmetric (prescription independent) formulation 
of this type of process. We show that, at the core of 
SUSY is the non-trivial definition of backward stochas- 
tic processes and the equilibrium distribution reached at 
long times. 

Multiplicative white-noise stochastic processes have 
a wide spectrum of applications, not only in physical 
and chemical systems [ll|, but also in biology [11] 
and even in evolution of economic variables [l3| . The 
paradigm of multiplicative stochastic dynamics is rep- 
resented by a Langevin equation with state-dependent 
noise, given in its simplest form by 



x{t) is a single stochastic variable and (^{t) is a Gaus- 
sian white noise, (C(t)) = 0, (C(i)C(i')) = ^{i - The 
drift force f{x) and the diffusion function g{x) are, in 
principle, arbitrary smooth functions of x{t). It is well 
known that the complete definition of the above stochas- 
tic differential equation is achieved only upon fixing an 
interpretation for the product g{x(t))C,{t) or, more tech- 
nically, upon the proper definition of the Wiener inte- 



gral There are several prescriptions used to de- 

fine this product, which can be summarized in the so 
called "generalized Stratonovich prescription" [T^ or "a- 
prescription" [l5| . This definition contains the usual Ito 
(a = 0) and Stratonovich (a 1/2) prescriptions as par- 
ticular cases. The interpretation of cq. (Il} depends on 
the physics behind a particular application. Once the 
interpretation is fixed, the stochastic dynamics is unam- 
biguously defined. 

The Fokker-Plank equation (FPE) associated with 
eq.(IT]) can be written as a continuity equation dtP{x, t) -\- 
dxJ{x,t) = 0, supplemented with the initial con- 
dition P{x,Q) — Pin(x), where the probability cur- 
rent is J{x,t) = [f{x) - (1 - a)g{x)dxg{x)]P{x,t) - 
{g'^{x)/2)dxP{x,t). We assume that the system con- 
verges, at long times, to an equilibrium probability dis- 
tribution, Pcq{x) — N exTp{—Ucq{x)) , computcd as a sta- 
tionary solution of the FPE with zero stationary cur- 
rent probability Jstix) = limt^oo Jix,t) — 0. The 
unique solution is [/oq (a;) = — 2 / {f{x')/g'^{x' j)dx'-\-{l — 
a)\ng'^{x). When considering / as a conservative force, 
obtained from a potential V{x) as / = —{l/2)g^dxV, the 
equilibrium distribution has the simpler form. 



C/eq(x) = U(a;) + (1 - a) Ing^x). 



(2) 



dx{t) 
dt 



f{xit))+gixit))at). 



(1) 



Therefore, the equilibrium distribution depends, not only 
on the independent functions f{x) and g{x), but also on 
the prescription used to define the stochastic process. 

Note that, only for a = 1 (Hanggi-Klimontovich pres- 
cription [H, [13]), we obtain the usual thermodynamic 
equilibrium distribution Ucq{x) = V{x). Provided we 
know that, for a given physical system, the equilibrium 
distribution is of the Boltzmann type with a previously 
known potential V{x), it is possible to determine the 
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stochastic process that converges to this distribution. To 
do this, the function f(x) is modified by adding a "spu- 
rious drift term", / — > / — (1 — ct)gdxg In this 
way, any correlation function is uniquely determined and 
equilibrium properties do not depend on the value of a 
(see, for instance, refs. [l^-[2l| for a discussion on the ex- 
perimental determination of f{x) and the spurious drift 
terms in nanoscopic systems). 

In this letter, we will adopt an alternative point of 
view. We consider a stochastic process modelled by 
eq.(IT]), with a particular value of < a < 1. In this con- 
text, the value of a is part of the model and the equilib- 
rium distribution reached at long times is given by eq. ([2|) . 
In this way, we can model stochastic dynamics whose fi- 
nal equilibrium distribution is not necessarily dictated h 
thermodynamics, such as econophysics applications 
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II. TIME REVERSAL AND EQUILIBRIUM 
PROPERTIES 

A key concept in our approach is the proper definition 
of time reversal evolution. The a-prescription is implic- 
itly associated with a particular direction in time. For 
instance, the Ito interpretation, also called the "prepoint 
discretization rule" , states that in a discretized time for- 
malism, the Wiener integral / g{x{t))£^{t)dt should be 
computed by evaluating g{x{tj)), using the value of x{t) 
at the "initial" time of each interval [tj,tj^i]. On the 
other hand, in the Hanggi-Klimontovich prescription, or 
"postpoint discretization rule", we evaluate (7(a;(tj+i)) 
taking the value of x{t) at the "final" time of each in- 
terval. Of course, in a backward evolution the role of 
"initial" and "final" are interchanged. In fact, the post- 
point prescription is sometimes called "interpretation in 
the backward sense" [13] ■ In the a-prescription, the dis- 
sipation function is evaluated taking a weighted aver- 
age of the "initial" and "final" values of x{t), g[x{Tj)] — 
g[{l — a)x{tj) + ax{tj^i)]. When considering a backward 
evolution (tj o tj+i), .g[a;(Tj)] is obtained by replacing 
a by 1 — a. Thus, if the forward stochastic process is 
performed in a definite prescription a, the backward tra- 
jectories evolve with the dual prescription 1 — a. Taking 
into account that the equilibrium distribution Ucq (given 
by eq.([21)) must be the same whenever it is reached in a 
forward {t — >• -l-oo) or backward {t ~oo) process, the 
correct definition of the time reversal transformation is. 
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Since 1 — 2q; is odd under the transformation a — >■ 1 — a, 
= /, as it should be. In the Stratonovich prescription 
{a = 1/2) , Tf = f and T simply corresponds to a change 
in the sign of the velocity. In any other prescription, the 
definition of the time reversal operator is more involved. 

It is interesting to see how the time reversal transfor- 
mation is implemented in the path integral formalism. 



We assign to each stochastic trajectory x{t), with end- 
points {xi,ti) and {xf,tf), a weight V{x{t)\xiti, x ft f) = 
det'\g) e^^I^l. The "action" S[x] is given by Q 



S 



dt 



2g2 



dx 



/ + agd^g 



(4) 



The conditional probability of the system being in state 
x/ at time t f , provided it was initially in state Xi at time 
ti, is given by P{xf,tf\xi,ti) = J Vx P{x{t)\xiti,Xftf) 
with the boundary conditions x{ti) = Xi and x{tf) = Xf. 

Using the transformation ([3]), we compute the 
time reversed conditional probability P{xi,tf\xf,ti) ~ 
TP{x f ,tf\xi,ti), which can be written in the path inte- 
gral formalism as P(xi, t/|x/, i^) = J Vx Aet~^{g) e~^^^\ 
with the boundary conditions x{ti) = Xf and x{tf) — Xi 
and the time reversed "action" S[x] — TS[x]. Notice 
that, in general, TP{xf,tf\xi,ti) ^ P{xi,tf\xf, ti), since 
the backward stochastic process {f,g,l — a) is not the 
same, but the dual of the forward one (/, g, a). 

After some tedious but direct manipulations, it is pos- 
sible to write the forward action S (eq.(j3])) as 5 = 
^ (y + a Ing^) I'^-fS*, in such a way that ,5 is time rever- 
sal invariant, S — TS. The time reversed action takes the 
form S ^ TS ^ -\ {y + {2 ~ Za) In g^)\\' -f S. Thus, 
the variation of the action under a time reversal trans- 
formation is just a total derivative term (for arbitrary 
values of a), S — S = Ueqixf) — Ueqixj), where Ucq{x) 
is the equilibrium potential of eq.© [23|. This property 
immediately implies the detailed balance relation, 

Pcq{Xi) P{Xf,tf\Xi,ti) = P{xf,ti\Xi,tf) Pcq{Xf), (5) 



where Pcq{x) 



Nexp {-Ucq{x)). Eq. ^ is valid for 
any prescription a and it reduces to the usual detailed 
balanced relation for additive noise when a = 1/2 since, 
in this case, P{xf,ti\xi,tf) — P{xf,ti\xi,tf). 
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In the context of stochastic thermodynamics 
can compute the increase of entropy in the medium 
associated with one specific trajectory as [1^ [26| . 

Asjrt[x{t)] = In P(x{t)\xiti,Xftf)/P(x{t)\xitf,Xfti)y 
For each stochastic trajectory beginning and ending in 
states sampled with the distribution p{x), the total en- 
tropy production is As = lnp{xi) — lnp(x/) — Ucq{xf) + 
Ucqixi). In the absence of an explicit time-dependent 
driving force, the stochastic entropy is a state function 
which depends only on the initial and final states. More- 
over, if we prepare the initial and final states with the 
equilibrium distribution p{x) = Pcq{x), we immediately 
conclude that As = 0, for each stochastic trajectory and 
for arbitrary values of a. The relevance of the "potential" 
Ucq in the entropy production was recently recognized in 
the context of the energy representation of a Brownian 
system in the Stratonovich prescription (27| . Here, we 
show that this concept is much more general, and it does 
not depend on the specific value of a. 
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III. SUPERSYMMETRIC FORMULATION 



sented as a linear transformation in the extended space, 



N-point correlation functions can be computed by us- 
ing the generating functional 



Z{J) ^ j det~^(g)e" 



s[x]+jr^dt'j{t')x{t') 



(6) 



J{t) is a source with compact domain, i. e. , it adia- 
batically goes to zero away from an interval (ti,tf) in 
which we will compute the correlation functions. The 
action S (eq.(j4l)) is computed over the entire time axis 
(—00 < t < +00). To deduce eq.®, we have sam- 
pled the initial state with the equilibrium distribution 
and we have assumed an ergodic stochastic evolution, 
Peq{xi) = \\m.T^oo P(xi,ti\x-T,—T). lu this context, 
total derivative terms in S[x\ do not contribute to the 
dynamics of any observable. In this way, a system de- 
scribed by equation ^ is automatically invariant under 
time reversal, since S and S — TS differ just in a to- 
tal time-derivative term. Moreover, if we impose the 
constraint a:(-|-oo) ~ x(— 00), the action is truly time- 
reversal-invariant, S — S — S. 

The complex structure of S, and the fact that, for gen- 
eral a, the usual rules of calculus do not apply, make 
calculations rather cumbersome. However, a simplifica- 
tion is possible by extending the space of trajectories, 
adding one auxiliary "bosonic" variable ip{t) and a cou- 
ple of Grassman variables ^{t), ^{t). The generating func- 
tional (eq.®) can be written in the extended functional 
space as 



Z[J]= / VxVipV^V^ e 



-S[x,cp,i,Q + J dtj{t)x{t) 



(7) 



where the action 5, in the new variables, is given by 

/' + 00 ( rl ^ 

I dti -m^^m + nmtm + ^^(trgi^r 



i(p{t) 



(8) 



Of course, through functional integration over the aux- 
iliary variables, we obtain the generating functional 
with the action given by Eq. Q. The a pa- 
rameter appears in the definition of the equal-time re- 
tarded Green's function (^(i)^(i))R = a. It should 
be noted that, for consistency, the equal-time advanced 
Green's function should be defined with the dual pre- 
scription {^{t)^{t))A = 1 - a. In this way, {^{t)£,{t))R - 
{£.{t)£,{t))A = 1 — 2a is odd under time reversal. The 
complexity of the stochastic calculus, associated with the 
definition of the Wiener integral, is now codified in the 
structure of the Grassmann variables. In this way, we 
can formally work out any calculation without explicitly 
indicating a specific prescription. 

The time- reversal transformation, eq.Q, is repre- 



r 




x{-t) 

vi-t) - f i{-t) 
fji-t) 

-vi-t) 



(9) 



where ry, fj are the time reversal Grassman variables. 
There is a subtlety in this transformation which involves 
boundary conditions. While \imt^±oo ^{t)S,{t) = a, the 
time reversed variables satisfy limi_>.-|-oo = 1 — a. 

Computing now the variation of the action under 
a time-reversal transformation, we obtain S — S = 
Ucq{xf) — Ucq{xi) + As, whcrc the entropy production 
associated with each trajectory (in the extended space) 
is 



As = 



(10) 



Interestingly, we observe that the entropy production 
vanishes due to fermionic number conservation, or more 
precisely, due to the invariance of the action under global 
phase transformation of the Grassman variables. 

The path-integral formalism is useful to make evident 
(otherwise hidden) symmetries of the stochastic process. 
For instance, the action given by eq.® is invariant under 
the transformation, Sx — A^, 5^ = 0, 6^ = iXip, 6ip = 0, 
where A is an anticommuting parameter. This nilpotent 
transformation {S^ — 0) is the famous BRS [1^ symme- 
try, discovered in the context of quantization of gauge 
theories. In the present context, it simply enforces prob- 
ability conservation, Z{0) — 1. There is another set of 
important symmetries related with equilibrium proper- 
ties that, together with BRS, is called supersymmetry. 
To display it explicitly, it is convenient to introduce a 
"natural" response variable (p = g'^ (^ip + id^ In g^CC) ^-nd 
to re-scale the Grassman variables 77 = g^, fj — g^, in such 
a way that the transformation {x,ip,^,£,) — > (x, (p,r],fj), 
has a trivial Jacobian, VcpVS^V^ — VipDrfDrj. In these 
variables, it is simpler to formally compute response func- 
tions. To see this, we slightly perturb the system out of 
equilibrium V —> V + h{t)x, and compute the linear re- 
sponse 



R{t,t') 



5{x(t)),, 



Sh{t') 



{x{t)iip{t')), (11) 



h=0 



where the last correlation function should be computed 
with the action S{x,(p,rj,ff). 

It is convenient to collect the transformed dynamical 
variables in the definition of a simple scalar superfield. 



$(t, e, 9) = x{t) + ei]{t) + f]{t)e + ip>{t)ee, 



(12) 



where we have introduced two "temporal" Grassman 
variables and 9. The system can be described in terms 
of two superpotentials V{^) and r(<I>), where the "diffu- 
sion potential" is introduced as g^^{x) = dx^{x). Then, 
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the time-reversal-invariant action S, can be re- written in 
terms of 9, 9) as 

J dtdOdO |L)r[$]£ir[$] + , (13) 

where we have defined the covariant derivatives D = dg, 
D (l/2)9g - 9dt, satisfying D"^ ^ D'^ = and 
{D,D]^~dt. 

The action of eq. (|T3| (and, of course, that one of 
eq.®) has a SUSY whose generators are 

Q = dg , Q = {ll2)de + 9dt , {Q, Q} ^ dt . (14) 

The graded algebra {Q,D} = {Q,D} = {Q,D} = 
{0, 5} = Q2 = = 0, guarantees that the variation of 
the action is a total derivative. 

Each of the three generators of SUSY {Q,Q,dt} im- 
poses several non-perturbative constraints on correlation 
functions, dt induces time translations and means that 
any two-point correlation function depends on time dif- 
ferences. On the other hand, Q, the generator of 9 trans- 
lations, represents the usual BRS symmetry, responsible 
for the probability conservation. Moreover, the invari- 
ance of the action generated by Q imposes, for instance, 

($(il,01,0l)$(t2,e2,^2)) = 

($(ii + e9i,9i + e/2, 9i)<P{t2 + e^2, ^^2 + e/2, ^2))(15) 

where e is an arbitrary Grassmann parameter. Eq. ljlSp 
implies, when written in components, the fluctuation dis- 
sipation relation R{t,t') - dt{x{t)x{t'))9{t - t'). 

For additive processes, the diffusive potential is linear, 
T{x) ^ X, and Eq.(fT3| reduces to the usual action de- 
fined with a single superpotential V{^). In this case, the 



tadpole theorem [9| guarantees that the stochastic evo- 
lution does not depend on a. However, multiplicative 
processes (non-linear F), induce derivative couplings in 
the superfield. These couplings are responsible for the 
a-dependent evolution that leads to the equilibrium dis- 
tribution of Eq.([2]). 

IV. SUMMARY AND CONCLUSIONS 

Summarizing, we have shown a hidden SUSY in a gen- 
eral multiplicative white-noise stochastic process. SUSY 
does not depend on the particular definition of the 
Wiener integral (a) and encodes equilibrium properties, 
such as fluctuation-dissipation relations. The main in- 
gredient in the SUSY formulation is the proper definition 
of time reversal and the a-dependent equilibrium distri- 
bution. Finally, we presented a covariant formulation 
based on two superpotentials related with the drift force 
and with the diffusion function. The SUSY structure is 
particularly useful to study multiplicative white-noise dy- 
namics, since any specific computation can be worked out 
in a unified way, without the specification of a particular 
value of a (nor a particular rule of calculus). Moreover, 
the nonperturbative constraints imposed by SUSY help 
to select the relevant diagrams that should be taken into 
account in a perturbative calculation. 
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